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ABSTRACT 

We present a new two-color algorithm, the "Stellar Bump Sequence" (SBS), that is optimized for 
robustly identifying candidate high-redshift galaxy clusters in combined wide-field optical and mid- 
infrared (MIR) data. The SBS algorithm is a fusion of the well-tested cluster red-sequence method of 
Gladders & Yee (2000) with the MIR 3.6/im - 4.5/im cluster detection method developed by Papovich 
(2008). As with the cluster red-sequence method, the SBS identifies candidate overdensities within 
3.6/im - 4.5/im color slices, which are the equivalent of a rest-frame 1.6/im stellar bump "red-sequence" . 
In addition to employing the MIR colors of galaxies, the SBS algorithm incorporates an optical/MIR 
(z' - 3.6/im) color cut. This cut effectively eliminates foreground 0.2 < z < 0.4 galaxies which have 
3.6/im - 4.5/im colors that are similarly red as z > 1.0 galaxies and add noise when searching for high- 
redshift galaxy overdensities. We demonstrate using the z ~ 1 GCLASS cluster sample that similar 
to the red sequence, the stellar bump sequence appears to be a ubiquitous feature of high-redshift 
clusters, and that within that sample the color of the stellar bump sequence increases monotonically 
with redshift and provides photometric redshifts accurate to Az — 0.05. We apply the SBS method 
in the XMM-LSS SWIRE field and show that it robustly recovers the majority of confirmed optical, 
MIR, and X-ray-selected clusters at z > 1.0 in that field. Lastly, we present confirmation of SpARCS 
J022427-032354 at z — 1.63, a new cluster detected with the method and confirmed with 12 high- 
confidence spectroscopic redshifts obtained using FORS2 on the VLT. We conclude with a discussion 
of future prospects for using the algorithm. 

Subject headings: galaxies: clusters and high-redshift - infrared: galaxies - cosmology: large-scale 
structure of universe 
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1. INTRODUCTION 

Massive galaxy clusters are important astrophysical 
objects. As the largest virialized structures in the uni- 
verse, their space density, mass function, and clustering 
provide rob ust constraints on cosm ological models (e.g. , 
Edge et alj Il990t I Eke et all 119981: iBahcall est all 120031: 
Gladders et all 120071: IVikhlinin et all 120091: iRozo et al l 
20101 : ISehgal et all 120111 : iBenson et all 120111 1 As the 
highest-density regions in the universe, they are also 
important cosmic laboratories for studying h ow envi- 
ronm e nt affects the evolu t ion of galaxies (e-g., iDressleil 
19801: iBldogh et alJ fl999l iPoggianti et all fl999t l200l 
Patel et al l l2009allb1: IVulcani et all 120101: IWetzel et all 
20111: IM uzzm et all 120121 and numerous others). For 
these reasons there has always been strong motivation 
to discover larger samples of clusters, increasingly 
more massive clusters, and naturally, clusters at the 
highest-possible redshifts. 

Substantial progress on all these fronts has been 
made in the last decade, particularly in finding clusters 
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to the highest redshifts. Ten years ago only a few 
clusters at z > 1 had been discovered and confirmed 
(jStanford et all 119971 : iRosati et all Il999h . Five years 
ago only a few clusters at z > 1.4 had been confirmed 
(jStanford et all 120051 . 12001 . Now, using a multitude of 
techniques, the number of confirmed clusters above z > 
1 is approaching ~ 100, and the number of confirmed 
clusters at z > 1.4 is now > 10. 

The progress in finding high-redshift clusters has 
occurred at an extremely rapid pace, but perhaps 
what has been most astonishing is that new high- 
redshift clusters are now being discovered routinely 
using different techniques and observations at a wide 
range of wavelengths. In previous decades searches for 
the highest-rcdshift clusters were almost completely 
dominated by X- ray surveys (e.g., iGioia et all Il990l : 
IRosati et al.lll998l ). whose superior sensitivity over wide 
areas allowed them to robustly discover massive clusters 
with relative ease. In recent years the advent of wide- 
field near- and mid-infrared imaging (hereafter NIR and 
MIR, respectively), as well as millimeter- wavelength 
surveys for clusters using the Sunyaev-Zel'dovich effect 
(hereafter SZE) has made these methods competitive 
with the latest generation of X-ray surveys. 

The current large sample of z > 1 c l usters 
comes from X-ray surveys (e.g. iMullis et all 12001 
Stanford et~all l2006t iBremer et all 120061: iPacaud et al l 



20071: iHenry et all 1201 Ot iFassbender et all l2011allbl: 



Nastasi et al. T201TI ISantos et all 1201 lh. space- 
based MIR surveys (e.g., IStanford et all 120051 
lEisenhardt et al.ll2008l : iMuzzin et all 2009; Wilson et all 
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2009t IPapovich et al.l l2010t iDemarco et all 120101: 
Brodwin et al.ll201lHStanford et al.ll2012tlZeimann et al. 
20121) ground-bas ed NIR surveys (e.g.. iGobat et al* 
201 it i lSpitler et al.ll2012D. millimet e r SZE surveys (e.g. 
Brodwin et al.l 120101 : iFolev et al.l 120111 : IStalder et all 
2012T ). as well as narrow-band and MIR searches around 
high-redshift radio galaxies (e.g.. IVenemans et al.l 120071 : 
iGalametz et"aT1l2010l |20H . 

Naturally, this array of detection methods has pro- 
duced a heterogeneous sample of high-redshift clusters. 
While this may create some difficulty when comparing 
the properties of clusters and cluster galaxies between 
various samples, there are also distinct advantages to 
having samples of clusters selected in diverse ways. 
For example, SZE and X-ray surveys detect clusters 
based on their baryonic mass in the form of hot gas, 
and therefore preferentially select halos with higher gas 
densities and higher gas fractions. Optical/IR surveys 
detect clusters based on their stellar mass, and therefore 
preferentially select halos with higher stellar-mass frac- 
tions. If the relative fraction of baryons in the form of 
gas and st ars varies substa ntially from halo to halo (e.g., 
iLin et al.l [20031 : iDai et aljfeoiof ). then these techniques 
are complementary - indeed almost necessary - for 
finding mass-complete samples of clusters. 

For technical reasons, X-ray and SZE cluster surveys 
have historically covered large volumes searching for 
the most massive systems, whereas optical/IR surveys 
have covered smaller volumes and have primarily discov- 
ered low-mass systems. Optical/IR surveys over large 
volumes are advantageous because they can discover 
complementary samples of the rare massive clusters 
typically selected in SZE and X-ray surveys. It is 
challenging to obtain optic al/IR data in multiple filters 
over wide areas; however, iGladders fc Yeel (l2000l 120051) 
demonstrated that clusters can be selected in wide-area 
surveys with a low contamination rate with just two 
photometric bands up to z = 1 using the red-sequence 
of cluster early-type galaxies as a high-contrast signpost 
for clusters. This red-sequence method was adapted 
into t he MIR for the SpARCS survey (|Muzzin et al.l 
120091: IWilson et all I2009D. whe re the 50 deg 2 Spitzer 
SWIRE fields (lLonsdale elaTl 120031) were observed in 
the z'-band and clusters to z ~ 1.4 were detected using 
a z' - 3.6/im color to identify galaxy overdensities. 

SpARCS was a logical extension of the red-sequence 
method into the MIR in order to search for clusters 
at higher redshift than optical surveys; however, the 
technique does require deep optical data to compliment 
the MI R data, a consid erable investment in observing 
time. IPapovichl ((2008) proposed a novel method for 
selecting high-redshift clusters using MIR data alone, 
thus avoiding t he ne ed for additional deep optical 
data. IPapovichl ()2008l ) demonstrated that high-redshift 
galaxies tend to have red 3.6/im - 4.5/im colors, and 
that because of the negative k-correction in these 
bands, even distant galaxies are bright enough to be 
detected in relatively shallow MIR imaging. He devised 
a color cut (3.6^m - 4.5^m > -0.1 AB) that allowed 
the identification of numerous cluster candidates in the 
SWIRE field at z > 1.3, and this approach led to the 
identificati on of the highest-re dshift known cluster at 
that time (jPapovich et al.ll2010D . 

In this paper we present the "Stellar Bump Sequence" 



(SBS) which is an improved algorithm for efficiently 
detecting high-redshift clusters in combined medium- 
deep optical/MIR imaging data. The method is a two 
color (three filter) algorithm based on combining the 
red-sequence cluster detection method used in SpARC S 
with the MIR method as developed by IPapovi ch (2008); 
IPapovich et al.l (|2010l) . In this paper we show that the 
hybrid of the two methods improves the reliability of 
candidate overdensities and also provides good-quality 
photometric redshifts for all cluster candidates. We 
apply this new technique to the XMM-LSS SWIRE 
field, and show that it recovers the majority of known 
clusters at z > 1 in that field, with good estimates of 
their photometric redshift. Lastly, we present a new, 
rich cluster spectroscopically confirmed at z = 1.63 that 
was discovered using the method. 

The layout of this paper is as follows. In §2 we present 
the SBS algorithm and demonstrate that the stellar 
bump sequence is clearly seen in spectroscopically- 
confirmed clusters at z ~ 1. In §3 we apply the 
algorithm to the XMM-LSS field and show that it recov- 
ers the majority of known clusters at z > 1 in that field. 
In §4 we present the confirmed cluster J022427-032354 
at z =1.63, a demonstration of the effectiveness of 
the algorithm for selecting high-redshift clusters. We 
conclude with a discussion in §5. Throughout this paper 
we assume a 17a = 0.7, Q m = 0.3, and H = 70 km 
s _1 Mpc" 1 cosmology. All magnitudes are in the AB 
system. 

2. THE STELLAR BUMP SEQUENCE ALGORITHM 

2.1. Motivation 

Photometric redshifts for galaxies are typically deter- 
mined by fitting the observed colors of galaxies with 
a range of templates for galaxy spectra. The features 
in galaxy spectral energy distributions (SEDs) that 
provide the tightest constraints on their photometric 
redshift are those that are the sharpest, particularly the 
Lyman break at 912A, the Balmer break at ~ 3800A, 
and the 4000A break in old stellar populations. An 
additional strong feature is the so-called "stellar bump" 
feature at rest- frame ~ 1.6pm which is caused by the 
minimum o pacity of the H _ ion at this wavelength (e.g., 
I John! 119881 ). This results in a maximum of the flux 
density at 1.6 um for galaxies whose ages are > 10 Myr 
(Sawicki 2002). As an illustration of the stellar bump 
feature, in Figure 1 we plot the SEDs of a single-burst 
simple stellar population (SSP) with solar metallicity 
formed at z = 10 as it would be seen at z = 0.5, 1.0, 
1.5, and 2.0. The S EDs h ave been generated using the 
iBruzual fc Charlotl (|2003l ) models. We note that we 
have used an SSP as a fiducial model for the purpose 
of illustration, but that any stellar population with an 
age > 10 My r also shows the 1.6/um feature (see e.g., 
ISawickill2002l IPapovichl 120081 ). The 1.6/zm stellar bump 
feature can be seen clearly in the SEDs in Figure 1, and 
its prominence resulted in it being discussed prior to the 
launch of Spitzer as a potentially powerful photometric 
redsh ift indicator ([Simpson fc EisenhardtJ Il999l : ISawickil 
12001) . 

In practice, the stellar bump does provide constraints 
on the photometric redshifts of galaxies; however, these 
tend to be much weaker than those obtained from the 
bluer features such as the 4000A break. In the lower 
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Fig. 1. — Left Panel: The evolving SED of an SSP with a Zf = 10 generated with the BC03 models as it would be observed at various 
redshifts. Filter response curves for the z', 3.6/im, and 4.5/im bands are shown below and the central wavelengths are marked with the 
dotted vertical lines. Top Right Panel: The z' - 3.6/im color as a function of redshift for the SEDs in the left panel. The z' - 3.6/im color 
increases monotonically with redshift and hence provides an estimate of the photo-z for an SSP. Bottom Right Panel: The 3.6/im - 4.5/im 
color as a function of redshift for the SEDs in the left panel. Due to the more complex SED shape in the rest-frame NIR, the 3.6/im - 
4.5/im color does not uniquely map onto redshift; however, if a z' - 3.6/im > 1.7 color cut is applied to cull low-redshift galaxies (see top 
right panel), the 3.6/tm - 4.5/im color becomes a viable estimate of the photometric redshift as those filters move through the rest-frame 
1.6/xm stellar bump feature. 



right panel of Figure 1 we plot the 3.6/im - 4.5/im color 
of this model as a function of redshift. The observed 
color maps uniquely onto redshift within the redshift 
interval 0.7 < z <1.7; however, the complete range in 
color is small (~ 0.7 mag), making the color only a 
modest-precision estimator of photometric redshift. In 
the upper right panel of Figure 1 we plot the observed z' 
- 3.6/im color of the model as a function of redshift. This 
color spans the 4000 A break at z > 1, and in contrast to 
the 3.6/im - 4.5/im color, increases monotonically with 
redshift at all redshifts, and has a substantial range 
(~ 4 mag). These properties make it a more powerful 
tool for identifying high-redshift cluster candidates and 
estimating their redshifts. 

In the SpARC S survey (jMuzzin et all 120091 : 
IWilson et al.l [2009) we used the z' - 3.6/im colors 
of galaxies to detect clusters up to z ~ 1.4. However, 
as Figure 1 shows, at z > 1.4 the z' - 3.6/im colors of 
the model become extremely red (z' - 3.6/im > 3.0). 
This implies that z'-band imaging that is substantially 
deeper than the MIR data are required to identify 
candidate high-redshift clusters with this method. This 
is approximately the limiting depth of the SpARCS z' 
data (z' < 24.2), which makes it challenging to select 
clusters at z > 1.4 robustly with the red-sequence 
method. 

IPapovichl (|2008f ) showed that provided only clusters 



at z > 1.3 were considered, the 3.6/im - 4.5/iin color of 
galaxies is sufficient to select high-redshift overdensities 
of galaxies. This novel idea allowed him to take full 
advantage of the deep SWIRE MIR data (a 5cr depth of 
~ 0.5 L* at z ~ 2), to select distant cluster candidates. 
IPapovichl (|2008l ) searched for overdensities of galaxies in 
position space that passed a single color cut of 3.6/im 
- 4.5/im > -0.1 in the SWIRE fields and discovered 
hundreds of z > 1.3 cluster candidate s including the 
confi rmed cluster C1G J0218.3-5010 (jPapovich et al.l 

This new method has been an important step for- 
ward for pushing MIR cluster detection to the highest 
redshifts, and C1G J0218. 3-5010 has been an important 
object for studying galaxy evolutio n in high-r e dshift , 
high-density environments (e.g.. iTran et al.1 120101: 
Papovich et all 120121: iQuadri et al.ll2012t iRudnick et al.l 



2012t iTadaki et al.ll2012D . Still, a single cluster almost 



certainly cannot be representative of the full population 
of clusters at z > 1.4 and much larger samples of clusters 
selected in the MIR that are mass-complete and have a 
high purity level would be extremely valuable. As of yet, 
it is unclear whe ther searching for clusters with only the 
IPapovichl (|2008f ) 3.6/im - 4.5/iin color cut is sufficiently 
robust to produce such samples. 

There are several disadvantages to the color-cut 
method as compared to red-sequence selection (e.g., 
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Fig. 2. — Left Panel: The observed 3.6/im - 4.5^tm color as a function of spectroscopic redshift for all galaxies in the FIREWORKS 
catalog of Wuyts et al (2008) with 3.6/im < 22.0. An SSP with a Zf = 10 generated using the BC03 models is overplotted for reference and 
is a reasonable description of the colors of most galaxies. This occurs because the 3.6/im - 4.5/tm colors of galaxies are fairly insensitive to 
their star formation history. Galaxies that lie considerably above the curve are likely to have a warm dust component heated by an AGN. 
Right Panel: Same as left panel, but with galaxies with z' - 3.6/im > 1.7 plotted. The color cut excludes most galaxies at z < 0.7 with red 
3.6/im - 4.5/tm colors and shows that for the remaining galaxies the 3.6/im - 4.5/im color can be used as a photo-z accurate to ± 0.15 for 
an individual galaxy. 



iMuzzin et all 120091 iWilson et all I2009D or multi-band 
photo metric redshift selection fe.g.. lEisenhardt et al.l 
l2008f ). Specifically, the fact that 3.6/im - 4.5/im does 
not uniquely map onto redshift means that galaxies at 
rcdshifts 0.2 < z < 0.4 have similar colors to those at 
z > 1.0 and add nois e when selecting high- redshift clus- 
ters. IPapovichl (|2008l ) argued that this contamination is 
minimized provided the MIR color cut was red enough 
(3.6/im - 4.5/im > -0.1); however, as we show in § 2.2.1, 
even with such a red color cut the contamination is 
considerable, particularly at bright magnitudes. 

The second issue is the lack of the use of additional 
color information beyond the initial 3.6/im - 4.5/im cut. 
Such selection compresses all structures at 1.3 < z < 
2.0 along the line-of-sight. This type of selection is 
similar to that adopted by the early-generation optical 
cluster surveys such as the Abell catalog. Those cluster 
catalogs suffer from low purity levels because numerous 
cluster candidates are in fact line-of-sight projections 
of low-mass groups. The projection effect problem in 
early optical cluster surveys has now been remedied by 
employing selection algorithms such as the red-sequence 
method, which detects clust ers in narrow color "slices" 
(jGladders fc Yed 120001 [20051) . hence reducing the prob- 
ability of lin e-of-sight p roject ions. Color slices are not 
used in the IPapovichl (|2008l) method, which suggests 
that a cluster sample selected using that method may 
suffer similar projection effect issues as in the early 
optical cluster surveys. 

The last issue with the MIR color-cut approach is 
that it does not provide information on the photometric 
redshift of cluster candidates (other than the cluster 
should be at z > 1.3). Whether they are used for 
cosmological purposes, or for studying the galaxy 
population, the photometric redshift distribution for a 
cluster sample is a key quantity. Furthermore, given the 
current lack of known clusters at 1.5 < z < 2.0, being 



able to differentiate cluster candidates at higher redshift 
compared to lower redshift is a distinct advantage for 
prioritizing followup observations. 

Here we outline some modifications that we have 
ma de to the MIR cluster detection method proposed 
by IPapovichl ()2008l ) that address some of these issues. 
The two main changes are the use of an additional opti- 
cal/MIR color cut, and the selection of clusters within 
multiple MIR color slic es. This ef f ective ly combines the 
MIR color approach of IPapovichl (|2008l ) with the well- 



tested clust er red sequence method of iGladders fc Yed 
(|2000l 120051 ). We show that the optical/MLR color cut 
substantially reduces the contamination of high-redshift 
galaxies by low-redshift galaxies and is so efficient at 
removing interlopers that the MIR selection method can 
be used down to much lower redshift {z ~ 0.85). Similar 
to the red-sequence approach, the 3.6/im - 4.5/im color 
slice method reduces the number of line-of-sight projec- 
tions, and also provides an estimate of the photometric 
redshift for each cluster candidate. 

2.2. Modifications to the Papovich Method 
2.2.1. z - 3.6/j,m color cut 

In the left panel of Figure 2 we plot the observed 
3.6/im - 4.5/im color vs. spectroscopic redshift for galax- 
ies with 3.6/im < 22.0 mag (the 5a depth of the SWIRE 
survey) in the GOO DS-South FIREWORKS catalog of 
iWuvts et al.1 (|2008l ). Overplotted as the solid line is 
the color of the single-burst model shown in Figure 1. 
The single-burst model is a reasonable description of the 
galaxy population at all redshifts, primarily because the 
3.6/im - 4.5/im colors of galaxies are insensitive to star 
formation history. 

Figure 2 shows that while there is a good correspon- 
dence between 3.6/im - 4.5/im color and redshift for 
galaxies at z > 0.7, there is a significant population of 
0.2 < z < 0.4 galaxies with colors similar to those at 



5 



SpARCS J003645-441050 SpARCS J02 1 524-034331 



SpARCS J161641+554513 SpARCS J022427-032354 



□ Spec Member 
Spec Interloper 





□ Spec Member 

□ Spec Interloper 



••• 




v.* 



z = 1.156 




z = 1.633 



17 18 19 20 21 22 17 
3.6jum 



19 20 



17 18 19 20 21 22 17 18 19 20 21 22 
3.6/xm 3.6/im 



Fig. 3. — The 3.6/im - 4.5/im vs. 3.6/tm color-magnitude diagrams for four SpARCS clusters. The three left-most clusters were discovered 
with the z' - 3.6/tm red-sequence method and have extensive spectroscopy from the GCLASS survey. Galaxies within R < 0.5 Mpc and 0.5 
Mpc < R < 1.0 Mpc are plotted as large and small symbols, respectively. Galaxies that have spectroscopic redshifts concurrent with being 
clusters members arc shown as red squares, those with spectroscopic redshifts inconsistent with being clusters members are shown as green 
squares. The cluster members belong to a fairly tight (scatter ~ 0.1 mag) sequence in color-magnitude space that becomes monotonically 
redder with increasing redshift. Similar to the red-sequence method, this "stellar-bump sequence" can be used as a signature to identify 
clusters in wide-field imaging. The solid lines and blue shaded regions are models for the stellar-bump sequence determined using BC03 
models. The 3.6/im - 4.5/rni CMD for SpARCS J022427-032354, which was identified from its stellar-bump sequence is shown in the 
right-most panel and is substantially redder than the lower-redshift SpARCS clusters. 



z > 1. iPapovichl (|2008l) argued that using a color cut of 
3.6/im - 4.5/im > -0.1 excludes most of the 0.2 < z < 0.4 
population. The number of z > 1.3 galaxies in the FIRE- 
WORKS spectroscopic sample with 3.6/im - 4.5/im > - 
0.1 and 3.6/im < 22.0 mag is 82. The number of z < 0.4 
galaxies matching the same color/magnitude criteria is 
23. Although the Papovich MIR color criteria does elim- 
inate some low redshift interlopers, clearly the interloper 
fraction is still no n- negligible, ~ 20%. We note that 20% 
is formally an upper limit for the contamination because 
the spectroscopy in GOODS-South is more complete at 
z < 0.4 than at z > 1.3. Still, the overall spectro scopic 
completeness is quite good (see lWuvts et al.ll2008l ). par- 
ticularly at these bright magnitudes, so the estimate of 
20% should be indicative of the actual interloper frac- 
tion. 

In the right panel of Figure 2 the FIREWORKS galax- 
ies are plotted again, this time excluding galaxies with 
z' - 3.6/im < 1.7. This removes the majority of the 0.2 
< z < 0.4 population without excluding the z > 1.3 pop- 
ulation. The color cut removes 20/23 of the low-redshift 
interlopers (~ 85%), but only removes 5 galaxies at z > 
1.3 (~ 6%). Therefore, with the additional z' - 3.6/im > 
1.7 criteria, the interloper fraction form galaxies at z < 
0.4 is improved by a factor of five to an impressively low 
~ 4%. 

Clearly a significant reduction in projection effects 
from low-redshift interlopers can be made by including 
a simple optical/MIR color cut in addi tion to the Pa - 
povich 3.6/im - 4.5/im color selection. IPapovichl (|2008l ) 
also showed similar reductions in the interloper fraction if 
a simple R-band magnitude cut was employed (R < 22.5 
mag). Indeed, making a color cut does not even require 
particularly deep optical data. The SpARCS z' data in 
the SWIRE fields is approximately 2.2 mag deeper than 
the 3.6/im data, deeper than is needed to make the z' - 
3.6/im > 1.7 cut. This implies that future surveys em- 
ploying the SBS selection could use ~ 3x less integration 
time than SpARCS in the z'-band, which would be only 
35 minutes per pixel on a 4m-class telescope. 



2.2.2. Color Slices vs. a Color Cut 

Within Figure 2, several overdensities in spectroscopic 
redshift are clearly visible. These overdensities have 
mean 3.6/im - 4.5/zm colors consistent with the single- 
burst model and fairly small dispersions around this color 
(~ 0.1 mag). This suggests that using slices in 3.6/im - 
4.5/im color would be an effective way to minimize pro- 
jection effects when selecting high-redshift clusters. This 
approach is identical to that taken for red-sequence se- 
lection, other than it uses very different rest-frame wave- 
lengths which bracket the 1.6/im stellar bump feature 
rather than the 4000A break. 

To further demonstrate the well-behaved nature of the 
MIR colors of cluster galaxies, in Figure 3 we plot the 
3.6/im - 4.5/im colors vs. 3.6/im flux for galaxies in the 
fields of three clusters at z = 0.869, 1.004, and 1.156 with 
extensive spectrosc opy from the GCLASS survey (see 
iMuzzin et al.1 12012). Confirmed spectroscopic members 
are shown as red squares and confirmed non-members 
are shown as green squares. Similar to the overdensities 
in FIREWORKS, the colors of the cluster members are 
nearly identical. It appears that similar to the red se- 
quence of early- types, there is a "stellar bump sequence" 
for galaxies in high-redshift clusters that increases mono- 
tonically with redshift (within this redshift window). 

The stellar bump sequence also appears to be quite 
tight. Overplotted in each panel is a model for the se- 
quence (see § 2.3) with a shaded region that represents a 
scatter of 0.1 mag. This is consistent with the rms scat- 
ter around the model for galaxies with 3.6/<m < 21.0 seen 
in the two lowest redshift clusters which are 0.093 and 
0.078 mag, respectively (rejecting >3er outliers). The 
data for those clusters is from the SWIRE survey that 
has 150s integration time per pixel and the median pho- 
tom etric error is 0.05 m ag. The z = 1.156 cluster (see 
also lDemarco et aLll2010T ) has data from a GTO followup 
program (PI: Fazio/ Wilson) and has an exposure time of 
2040s per pixel. The scatter in the stellar bump sequence 
for that cluster is substantially smaller, 0.047 mag (again 
rejecting >3er outliers). The median photometric error 
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in the deeper data is 0.02 mag. Subtracting the photo- 
metric scatters from the measured scatters in quadrature 
for all three clusters suggests that the stellar bump se- 
quence has an impressively small intrinsic scatter of 0.04 

- 0.07 mag. 

Figure 3 shows that selecting clusters in narrow slices 
of 3.6/im - 4.5^m color should be an effective way to min- 
imize projection effects. The sequences in Figure 3 are 
spaced such that there should be little overlap within the 
photometric scatter for clusters Sz > 0.15 apart in red- 
shift space. This is not as good as the precision from red- 
sequence selection where sequences typically only overlap 
within the photometric scatter for clusters with Sz < 0.05 
(jGladders fc Yedl2005t IMuzzin et al.ll2008ft : however it is 
a significant improvement over no color selection what- 
soever. 

Lastly, we note that the sol id lines shown in Figure 3 
are not fits to the data, but are lBruzual fc Chariot! (2003) 
models for the 3.6/mi - 4.5/im sequence at the redshift 
of the cluster assuming a single burst with a zt = 10 . 
These are remarkably good descriptions of the data and 
show that the color slice method should also return a 
reasonable estimate of the photometric redshift for the 
cluster (we discuss this further in § 5). 

2.3. The SBS Algorithm 

The SBS algorithm is nearly identical to the 
red-sequence d etecti on algorithm proposed by 
IGladders fe Yeel (12000ft a nd applied in IGladders fc Yeel 
(|2005ft and IMuzzin etaLl ([2008:) and we refer to those 
papers for a detailed step-by-step description of the 
method. In brief, the method is the following. Firstly, 
several SSP models for the apparent 3.6^m - 4.5/mi 
colors of g alaxies as a function of red shift is determined 
using the iBruzual fc Chariot! (|2003t ) models. These 
models are normalized using the U - V, V - I and J 

- K colors of the red sequence in the Coma cluster 
(|Bower et al.l 11992ft . These models are then used to 
construct color slices analagous to those plotted in 
Figure 3. For each color slice the probability that a 
given galaxy in the survey area belongs to that slice is 
calculated, the result being a probability weight for each 
galaxy, in each slice. In this process we only calculated 
weights for galaxies that pass the z' - 3.6/^m > 1.7 
color-cut. 

Probability maps are then created for each slice by 
placing galaxies on a pixelated grid of the survey field in 
their observed positions, weighted by their probabilities. 
The maps are then smoothed with a kernel with a 
physical size of 0.5 Mpc FWHM, and object detection 
is performed to identify sources above the background. 
Sources in the probability map are then flagged as 
candidate galaxy clusters. Sources are typically found 
in several color slices and a final catalog of candidates 
is made by connecting detections between slices and 
assigning candidates a photometric redshift based on 
the slice where they have the most significant detection. 
An example of SBS probability maps in the XMM-LSS 
SWIRE field in two different color slices is shown in 
Figure 4. 

3. APPLICATION TO THE XMM-LSS SWIRE FIELD 

The XMM-LSS field contains X-ray data from the 
XMM - Newton satellite taken as part of the XMM- 



LSS cluster survey ()Pierre et al.l I2OO60 . The field has 
been searched extensively for high-redshift clusters us- 
ing a variety of techniques, making it an ideal testbed 
for the SBS algorithm. There are four spectroscopically- 
confirmcd z > 1 X-ray-selected clusters from the XMM- 
LSS: XLSSC005, XL S SC029, XLSSC048 and XLSSC046 
(iBremer et al.l 120061: iPacaud et all 120071: lAdami et ail 
I2011J). There is one confirmed cluster from the SpARCS 
survey (SpARCS J021524-034331) discovered w i th th e 
z' - 3.6//m red-sequence method (|Muzzin "eTall[20Tl . 
The field also contains the spectroscopically-confirmed 
cluster CIG J0218.3 -510 at z = 1.62 discovered by 
iPapovich et al.l (|2010l ) using the MIR color-cut method. 

The field contains four channel IRAC im aging from the 
SWIRE survey (see lLonsdale et~ai1l2003l ) covering a to- 
tal area of 9.4 deg 2 . There is also deep z' data taken as 
part of the CFHTLS-Wide survey, as well as an extended 
region of z' data coverage taken as part of a CFHT fol- 
lowup program (PI: Pierre, ID: 2006BF11). We use these 
datasets to apply the SBS method and search for clusters 
at 0.85 < z < 2.0 

In Figure 4 we show the probability maps generated us- 
ing the SBS method. The top panel shows the color slice 
corresponding to a z p hoto = 1-0, the bottom panel shows 
the same field but generated with galaxies in the color 
slice Zphoto = 1-6. In both panels the width of the slice 
in color corresponds to a width in redshift of Sz ~ ± 0.1. 
In the z p hoto = 1-0 slice there are several overdensities 
that are clearly visible. Those that are associated with 
the confirmed z > 1 clusters are circled in red. When 
object detection is run, the SBS algorithm recovers 3/4 
X-ray selected clusters, as well as the red-sequence se- 
lected SpARCS cluster. There appears to be no overden- 
sity associated with the 4th X-ray cluster (XLSSC005), 
although we note that for this cluster the spectroscopic 
redshift is based on only 2 members an d therefore its 
spect roscopic redshift may not be secure (jPacaud et al.1 
12007ft . 

Figure 4 shows that the SBS selection appears to work 
well, recovering most of the high-confidence clusters at 
z ~ 1 selected with other methods. As shown by the blue 
circles in Figure 4, the method has also identified two 
other 2^1 overdensities with even higher significance 
than any of the confirmed clusters. These two candi- 
dates, SpARCS J022158-034000 and SpARCS J022056- 
033348 have estimated photometric redshifts of 1.08 and 
1.10, res pectively and are not in the latest XM M-LSS 
catalogs (jPacaud et al.l 120071 : lAdami et al.ll2011h . 

We note that none of the spectros copically-confirmed 
XMM-LSS clusters at 0.2 < z < 0.4 (lAdami et al.ll2011ft 
are associated with overdensities in the probability map, 
even though they have similar 3.6/im - 4.5/mi colors (see 
§ 5). This further highlights the effectiveness of the z' 
- 3.6/^m color cut in eliminating low-redshift interlopers. 
The removal of these interlopers also means that the SBS 
method can be robustly applied to redshifts as low as z ~ 
0.8, cr eating a larger r edshift range than can be probed 
by the iPapovichl (|2008ft color cut method. 

4. THE CLUSTER SPARCS J022427-032354 AT z = 1.6 

The bottom panel of Figure 4 shows the color slice that 
is centered at a z p hoto — 1-6. There appear to be fewer 
cluster candidates in the z = 1.6 slice compared to the 
z = 1 slice. This is partially because the z = 1.6 map 
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Fig. 4.— Top Panel: Smoothed density map of galaxies in the XMM-LSS SWIRE field with z' - 3.6^im > 1.7 and -0.5 < 3.6/xm - 4.5/im 
< -0.3. The field of view on the sky is ~ 3 X 2 degrees, north is up, east is to the right. The locations of spectroscopically-confirmed 
clusters discovered with X-ray observations (XLSSC clusters) and the red-sequence method (SpARCS J021524-034331) are shown with red 
circles. Those clusters are clearly associated with overdensities in the map and this shows that the 2-color SBS-method recovers clusters 
discovered with other search techniques. Shown in blue are two new cluster candidates discovered with the technique (see text). Bottom 
Panel: Same as top panel but for galaxies in a redder color slice, 0.0 < 3.6^tm - A.5fim < 0.2. The cluster SpARCS J022427-032355 is 
clearly visible in the upper right and is the most significant overdensity in the map. The cluster discovered by Papovich et al. (2010) using 
the SWIRE observations is also visible in the bottom left and is also one of the more significant overdensities in this color slice. 



is noisier than the z = 1 map due to the fixed depth of 
the observations. Even so, the data are deep enough (M 
< M* + 1.3) to robustly detect massive clusters so the 
reduced number of candidates suggests a significant de- 
cline in the number of massive clusters between z — 1.0 
and z = 1.6, consistent with exp ectations from ACDM 
models (e.g.. lHaiman et aT]|2001|) . 

The most prominent structure in the z = 1.6 slice is 
the candidate cluster in the northeast corner of the im- 
age which we designate SpARCS J022427-032354. Also 
visible in the southwest corner of the field is CIG J0218.3- 



510 discovered bv lPapovich et all (|2010t ). That cluster is 
also one of the more prominent overdensities in the color 
slice. 

In the left panel of Figure 5 we show a gz'3.6/xm 
color-composite of SpARCS J022427-032354. The clus- 
ter candidate is clearly visible as an overdensity of red 
galaxies in the center of the i mage. Similar to many 
other high-redshift clusters (e.g.. lEisenhardt et al.ll2008t 
iPapovich etal] 120101: iMuzzin et al.ll2012f ). the structure 
appears to be filamentary without a dominant central 
galaxy. 





Fig. 5.— Left Panel: gz'3.6/im color composite of the cluster SpARCS J022427-032354. The field-of-view is 7' X 5'. Right Panel: Same 
as left panel but with spectroscopically-confirmed cluster members shown in white and spectroscopically-confirmed foreground/background 
galaxies shown in green. 
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Fig. 6. — Histogram of spectroscopic redshifts in the field of 
SpARCS J022427-032354. The inset shows the distribution around 
the cluster redshift. 

In the right panel of Figure 3, we have also plotted the 
3.6/im - 4.5/im colors of galaxies near SpARCS J022427- 
032354 as a function of 3.6/im. Similar to the GCLASS 
clusters, SpARCS J022427-0323354 has a well-defined 
stellar bump sequence; however, that sequence is sub- 
stantially redder, consistent with a z p hoto = 1.6. Given 
its high significance and estimated redshift, SpARCS 
J022427-032354 was designated as a high-priority tar- 
get for spectroscopic followup in our vl.O SBS-selected 
cluster catalog of the XMM-LSS field. The cluster was 
observed using two multi-slit masks with the FORS2 in- 
strument on the VLT as part of Program 085.A-0613(B) 
(PI: Demarco). The on-sky integration time for each 
mask was 3.75 hrs. The 2D spectra were reduced using 
standard techniques for flat-fielding, wavelength calibra- 
tion, and rel ative flux norma l izatio n using the pipeline 
described in iDemarco etldl (|2005l 120071) . Once ID 
wavelength-calibrated spectra were extracted, spectro- 
scopic redshifts were determined both by cross correla- 
tion and eye-examination of faint spectra for emission 
lines. Full details of the procedure will be discussed in a 
future paper (G. Wilson in preparation). 

From the two spectroscopic masks a total of 56 high- 
confidence redshifts were determined. In Figure 6 we 
plot the spectroscopic redshift distribution of galaxies 
within 3 arcmin of SpARCS J022427-032354. There is 
a clear overdensity of galaxies at z = 1.63, with a total 
of 12 galaxies in the range 1.62 < z < 1.64 (see inset 
of Figure 6). In Figure 7 several examples of spectra in 
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Fig. 7. — Example spectra of cluster members of SpARCS 
J022427-032354. Portions of the spectra with strong atmospheric 
absorption have been masked out. The locations of the Mgll and 
Mgl absorption lines as well as the [Oil] emission line are marked. 
The majority of the spectroscopic members show strong [Oil] emis- 
sion and Mgll absorption suggesting active star formation. Shown 
is blue is the mean stacked spectrum of cluster members in a lower 
redshift cluster that is dominated by galaxies whose star formation 
is quenched. The top two spectra show much weaker [Oil] emis- 
sion and have a continuum shape that matches the lower-redshift 
cluster suggesting they may be quenched galaxies. 

this redshift range are shown. The majority of galax- 
ies in this redshift range with spectra show clear [Oil] 
emission (9/12); however, there are three spectra that 
have weak or possibly no [Oil] emission, two of which 
are shown in Figure 7. For these galaxies there is a 
nearby night sky line that makes it difficult to put strong 
constraints on the level of [Oil] emission; however, the 
continuum shape is similar to that expected in galaxies 
with quenched star formation. Overplotted in blue in 
Figure 7 is the mean stacked spectr a of galaxies in the 
cluster SpARCS J0 03550-431224 fsee I Wilson etlHl2059t 
iMuzzin et al.l [20121 ) which primarily contains quenched 
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Fig. 8. — Left Panel: The observed z' - 3.6/^m color of the red-sequence as a function of spectroscopic redshift for clusters within the 42 
square degree SpARCS survey. Points are color-coded based on the source of the spectroscopic redshift. Two SSP models from the BC03 
code with formation redshifts of 4.0 and 10.0 are shown and are good descriptions of the data. Right Panel: The observed 3.6/im - 4.5/tm 
color of the stellar-bump sequence as a function of redshift for the same clusters, with the same BC03 models. Both the z' - 3.6/tm and 
3.6/jm - 4.5/tm color-sequences of clusters can be used as an accurate photometric redshift estimate up to z ~ 1.7. 



galaxies. The UV continuum of the template matches 
the UV continuum of the weak [Oil] galaxies in Figure 
7 quite well, suggesting that these two galaxies may be 
quenched cluster members at z = 1.63. 

In the right panel of Figure 5 the positions of these 12 
cluster galaxies are plotted as white squares and the po- 
sitions of the remaining spectroscopic redshifts as green 
squares. Most of the galaxies with redshifts 1.62 < z < 
1.64 are concentrated in the central few arcmin of the 
field suggesting this is a cluster of galaxies as compared 
to a large, unvirialized "sheet" of galaxies. There are 
however several galaxies at the same redshift that ex- 
tend to ~ 1 Mpc, suggesting that this may be a massive 
cluster in the process of forming and hence the system 
may be only partially virialized. A more detailed analy- 
sis of the dynamics of the system will be presented in a 
future paper (G. Wilson, in preparation). 

We note that there are many red galaxies near the cen- 
ter of the galaxy distribution that do not have redshifts. 
Several of these were targeted in the spectroscopic masks 
but we were unable to determine a spectroscopic redshift 
due to the target being too faint in the observed optical. 
These galaxies have very red z' - 3.6/^m colors and could 
potentially be additional quiescent galaxies in the clus- 
ter core. If they are also at z = 1.633, then future NIR 
spectroscopy would be valuable to determine redshifts 
using the absorption features that are typically present 
in galaxies with old stellar populations. 

Overall, with 12 confirmed members, it appears that 
SpARCS J022427-032354 is the first confirmed cluster 
discovered using the SBS selection technique. 

5. DISCUSSION AND CONCLUSIONS 

The confirmation of SpARCS J022427-032354 shows 
that the SBS method is capable of selecting high-redshift 
clusters in combined optical/MIR datasets. Here we in- 
vestigate how the SBS photometric redshifts compare to 
those determined from the typical red-sequence photo- 
metric redshifts. In the left panel of Figure 8 we plot 



the median z' - 3.6/im colors of red-sequence galaxies 
from several different samples of clusters in the SpARCS 
fields as a function of spectroscopic redshift. The lowest- 
reds hift clusters come f rom the XMM-LSS cluster sam- 
ple ()Adami et al.l 1201 It ) and are X-ray selected. The 
GCLASS clusters are z' - 3.6/xm red-sequence se l ected . 
The "other" clusters come from iDemarco etall (l2010h 
(red-sequence selection) and lHashimoto et al.l ( 2005t ) fX- 
ray selection). Also shown are CIG J0218. 3-510 and 
SpARCS J022427-032354. Overplotted are two single- 
burst models with formation redshifts of Zf = 4.0 and 
Zf = 10.0. The models are good descriptions of the red- 
shift evolution of the color of the observed red-sequence 
and the rms scatter of the data compared to the model 
is 0.18 mag, implying that the z' - 3.6/im red-sequence 
color provides photometric redshifts for clusters to an ac- 
curacy of Sz — 0.08 from 0.1 < z < 1.6. 

In the right panel of Figure 8 we plot the median 3.6/xm 
- 4.5/im colors of the stellar bump sequence for the same 
clusters, as well as the same single-burst models. The 
model is also a good description of those data. The rms 
scatter in color is 0.05, which implies a redshift accuracy 
of Sz = 0.05 for clusters at 0.8 < z < 1.6. This is better 
than the z' - 3.6/im colors as a photometric redshift and is 
surprising given that 3.6/nm - 4.5/im spans a much smaller 
range over this redshift. The slightly better z p hoto may 
be because it is less sensitive to star-formation history, or 
because the MIR zeropoints are more stable than ground- 
based zeropoints over wide fields. Regardless, Figure 6 
demonstrates just how well the models describe the color 
evolution of the population and that good quality z p h t 
can be determined using the 3.6/^m - 4.5/xm colors of 
galaxies in candidate clusters. 

One caveat with the SBS method is that it is likely 
to have a higher false-positive detection rate compared 
to other optical cluster detection techniques such as the 
red-sequence method. Although the method does use 
3.6/im - 4.5/im color slices, the observed 3.6/xm - 4.5/xm 
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colors of galaxies do not evolve as rapidly with redshift 
as observed colors that span the 4000A break such as 
z' - 3.6/im. Therefore, line-of-sight projection will oc- 
cur more frequently. The false -positive rate for red- 
sequence selection i s ~ 5% (see, iGladders fc Yed 120051 : 
iGilbank et al.ll2007t ). Given the limited color range, the 
SBS false-positive rate may be as much as several times 
higher than this. At present, with only one confirmed 
cluster it is not possible to quantify the purity level of 
cluster samples selected with SBS. In the future, quan- 
tifying the false-positive rate will be important for using 
large samples of these clusters. 

Overall, the SBS algorithm appears to be an effec- 
tive means of identifying candidate high-redshift clusters 
in combined optical/MIR data. It is b ased on combin- 
ing th e cluster red-sequence met hod of [G ladders fc Yee 
(2000) with the MIR method of iPapovichl (|2008l ); how^ 
ever, the combination of the two provides distinct ad- 
vantages over each individually. Using the MIR colors of 
galaxies allows the cluster searches to continue to higher 
redshift over wider areas than is feasible using ground- 
based optical imaging. Likewise, the use of MIR color 
slices and the optical/MIR color cut reduces the low- 
redshift interloper fraction and the probability of line-of- 
sight projections which increases the purity level of the 
cluster sample. As shown in Figure 8, it also allows for 



an estimate of the photometric redshift for all clusters in 
the survey with good accuracy. 

The SBS algorithm has potential to provide a large, 
homogenously-selected sample of clusters in the coming 
years. We have already applied it to the 50 degree 2 
SWIRE fields and discovered several hundred cluster can- 
didates at z > 1 . In this paper we presented confirmation 
of the cluster SpARCS J022427-032354 at z = 1.63; how- 
ever, several other clusters discovered with the method 
in the SWIRE area have also been confirmed and will be 
discussed in future papers (G. Wilson, in preparation, 
R. Demarco, in preparation). In the future the method 
could also be applied to the available wide-field imaging 
datasets from Spitzer. At present, these cover of or- 
der several hundred square degrees and could provide an 
intermediate- mass cluster sample up to z ~ 1.7, some- 
thing which has potential for providing both interesting 
constraints on the growth of massive structures as well 
as the evolution of the cluster galaxy population. 
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